Exploration of a modified density dependence in the Skyrme functional 
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A variant of the basic Skyrme-Hartree-Fock (SHF) functional is considered dealing with a new 
form of density dependence. It employs only integer powers and thus will allow a more sound basis 
for projection schemes (particle number, angular momentum). We optimize the new functional with 
exactly the same adjustment strategy as used in an earlier study with a standard Skyrme functional. 
This allows direct comparisons of the performance of the new functional relative to the standard one. 
We discuss various observables: bulk properties of finite nuclei, nuclear matter, giant resonances, 
super-heavy elements, and energy systematics. The new functional performs at least as well as the 
standard one, but offers a wider range of applicability (e.g. for projection) and more flexibility in 
the regime of high densities. 
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I. INTRODUCTION 

Self-consistent mean-field models have grown to be a 
major tool in the theoretical description of nuclei and 
nuclear dynamics. There are several approaches to nu- 
clear SCMF from which the most widely used are the 
Skyrme-Hartree-Fock (SHF) model [3], the Gogny force 
0, and the relativistic mean-field model (RMF) 0, [|. 
All these models employ effective interactions which pro- 
vide reliable nuclear structure properties and low-energy 
excitations at the level of a mean-field description. The 
general strategy for constructing nuclear energy function- 
als for high-quality calculations is to deduce the formal 
structure from principle considerations Q and to rely 
on a phenomenological adjustment of the model param- 
eters. There is a long history of SHF development and 
optimization (for recent reviews see [H, @|) and yet, there 
is still demand for improvement as can be seen from the 
many recent updates of the parameterizations, for re- 
cent examples from SHF see e.g. 0-fHl- In spite of 
all these activities, there remain several open problems 
which touch the structure of present days SHF function- 
als fl2j . One such point is the way in which the density 
dependence is modeled. It used to start from adding to 
the zero-range two-body term a "three-body term" as- 
sociated with a density dependence oc p 3 where p is the 
local nucleon density [13|. This was found to be too re- 
strictive and has been generalized to a profile cx p 2+a 
with a ~ 0.16-0.3 in order to allow a better description 
of fission barriers and monopole oscillations ^4|. Since 
then, this form for the one density dependent term was 
carried forth unchanged. The reason is that most empir- 
ical information on nuclear structure is associated with 
densities near bulk equilibrium density po « 0.16 fm -3 . 



Systematic information in a wide range of densities has 
not yet been acquired contrary, e.g., to the case of model- 
ing liquid 3 He [Hill- 

Nonetheless, the traditional form 
oc p 2+a is an arbitrary guess, it is rather rigid what ex- 
trapolation to the low density region (=weak binding) is 
concerned, and in particular, the non-analytical behavior 
for non-integer a turns out to be a great hindrance for 
particle-number projection [12, Q3 • It is the aim of this 
paper to present and discuss an alternative form of the 
density dependence for the SHF functional which is built 
in the form of rational approximants using only integer 
powers of p. This will then provide a SHF functional 
which can safely be used in projection schemes, partic- 
ularly for particle- number projection [jjjj . Moreover, we 
will exploit the generalized form to explore the effect of 
varied density dependence on bulk properties of finite 
nuclei. The new SHF functional is calibrated and inves- 
tigated with the techniques of least-squares fitting and 
the basic data set as developed in a previous large scale 
study 0. This allows a direct comparison with the for- 
mer results using the traditional density dependence and 
it provides the same systematics of variations. 



II. THE ENERGY FUNCTIONAL 

The SHF functional is formulated in terms of the fol- 
lowing local densities and currents (for details see 0, lioj 
and appendix [Al : density p q , kinetic-energy density T q , 
spin-orbit density 3 q , current j 9 , spin density u q , and 
pair density £ 9 , where the index q £ {p, n} labels protons 
and neutrons; total densities are denoted as p = p p + p n 
and similarly for the other densities. Starting point is an 
expression of the total energy 
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From this total energy it is straight forward to derive the 
corresponding mean field equations by variation with re- 
spect to the single-nucleon wavefunctions and the pairing 
amplitudes. 



A. The Skyrme energy functional 

The Skyrme functional in the narrower sense is given 
by the Skyrme energy density £gk comprised in the terms 
(| lcil ITfj) . The parts (|ld[) and fTfJ carry the zero-range 
two-body terms including effective mass and spin-orbit 
couplings. All density dependence is collected in the 
term fTc]) . Its first contribution oc 63,63 represents the 
traditional density dependence. The second contribution 
adds the essential new piece where the density depen- 
dence is modeled as a rational approximant. Note that 
this term employs only integer powers of p. This new 
form has thus the great advantage that it complies with 
particle-number projection. The idea is that this rational 
density dependence replaces the traditional dependence 
p a . Thus we can switch the models by the coefficients. 
The traditional model is recovered by setting 63, b' 3 = 0. 
The new model with rational density dependence uses 
non-zero 63, 6 3 and switches off the traditional form with 
63,63 = 0. Both forms associate the same density de- 
pendence to isoscalar and isovector terms. One may in- 
troduce more isovector freedom in the density dependent 
terms by using all terms, i.e., 63,63 7^ together with 



63, 63 7^ 0, but fixing a = 1 to maintain the crucial fea- 
ture that only integer powers of p are used. We will 
consider this extended option to explore the impact of 
(isovector) density dependence. 

One may wonder why one needs such an involved ratio- 
nal form. A straight-forward approach would start with 
a mere Taylor expansion for the density dependence, e.g. 
oc p 2 [cip 1 + c 2 p 2 + csp 3 ...] . It turns out that this sim- 
ple form is not able to reproduce the appropriate density 
dependence at low densities. The form p a with small a 
performs, in fact, very well in this regime which explains 
its great success in the development of the SHF func- 
tional. It is clear that the non-analytical p a resists a 
straightforward Taylor expansion. The rational approxi- 
mant is much more flexible and allows a good fine tuning 
also in the low-density regime. 

The spin-orbit part of the SHF functional (fTfj) contains 
the tensor term with a switch factor ?y t i s . The importance 
of and the need for this term is still a matter of debate 
P. ITol. [l9| . The new parameterization developed here will 
employ the tensor terms, i.e. we deal with jftis = 1 • 

The Skyrme energy as given above contains the min- 
imum number of time-odd terms which are needed to 
make the functional Galilean invariant [2(|. There are 
many more time odd terms conceivable at the level of 
given functional P, [2l], [22| . They are unimportant for 
even-even nuclei and will be ignored here. 

The parameters 6i, b[ used in the above definition of the 
Skyrme functional are chosen to give a most convenient 
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formulation of the energy functional, the corresponding 
mean-field Hamiltonian and residual interaction. Tradi- 
tional expressions for the Skyrme functional use a dif- 
ferent labeling. The relation to the standard Skyrme 
parameters U, x\ 

HL H IM HI is S iven in appendix [Bj 

B. Additional terms 

The direct Coulomb term employs the mere proton 
density p p and not the charge density which would be 
obtained from folding with the nucleonic charge distri- 
butions The difference is, in principle, sizeable for 
heavy nuclei, but easily compensated for by the other 
terms in the functional. We keep the above form for rea- 
sons of simplicity. 

The Coulomb-exchange functional, the second term in 
( flgl , has been made switchable by the parameter f?ccx- 
Unfortunately, the LDA for Coulomb exchange invokes 
a non-integer power of density which, again, inhibits 
particle-number projection. Thus we will omit this con- 
tribution by using f?ccx = in the new parameterizations. 
Coulomb exchange is a small contribution and the miss- 
ing term will easily be taken up by the other terms in the 
functional. 

For the center-of-mass correction E cm , we use the stan- 
dard recipe as given in eq. (JTTJ) . Note the comment "after 
variation" therein. The E cra is, in fact, the expectation 
value of a two-body operator which makes the mean-field 
equations in fully variational treatment rather cumber- 
some. As E cm is a small correction, we evaluate it (non- 
variational) at the end of the mean-field calculation for 
then given single-particle wavcf unctions. 

The pairing functional (jlh[) contains a continuous 
switch, the parameter /?o,pair where volume pairing is 
recovered for po,pair — ► oo and surface pairing for 
Po, pah- = 0.16 fm~ 3 . We will use po,pair as a free param- 
eter and usually obtain a form of the pairing functional 
which stays in between the extremes of volume and sur- 
face pairing fTo| . 

C. Adjustment strategy and the variants in this 
survey 

The parameters of the new SHF functional are deter- 
mined by phenomenological adjustment as it is done for 
most SHF functionals. We employ exactly the same fit 
strategy and data as used in a recent large scale survey 
[Icij . The data set has been selected to cover spheri- 
cal nuclei with negligible correlation effects [25j which 
amounts to isotopic and isotonic chains of semi-magic 
nuclei (leaving out the mid-shell regions). As leading ob- 
servables we consider binding energy and key pattern of 
the nuclear charge distribution: r.m.s. radius, diffraction 
radius and surface thickness {26j . For pairing properties 
we include the three-point gaps as deduced from the odd- 
even staggering of binding energies and for determining 
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TABLE I: Definition of the parameterizations. The acronym 
RD stands for "Rational Density dependence" and RM for 
"Rational and More density dependence" . The entry "std" in 
the rightmost column indicates that the standard pool of fit 
data from finite nuclei is used fiol ) - Further entries, if present, 
indicate the nuclear matter properties which had been added 
to the standard set: a 3yrn = symmetry energy at bulk equilib- 
rium (fixed at 30 MeV), k = isovector sum-rule enhancement 
factor (fixed at 0.4), m* /m — effective nucleon mass (fixed at 
0.9). 



the spin-orbit parameters we take into account selected 
spin-orbit splittings in doubly-magic nuclei. The ac- 
tual adjustment is performed with straightforward least- 
squares techniques. These allow not only to determine 
the optimal parameters but also provide a measure for 
the uncertainty on the parameters and deduced observ- 
ables [13 [13. 

Table U shows the variants of the model studied in the 
following. The parameterizations SV-bas and SV-min 
rely on the standard SHF functional (note the 63, b' 3 = 0) 
and were presented in great detail in flpj ] . SV-min re- 
sults from a straightforward fit to the data from finite 
nuclei. The "min" stands for minimization of the qual- 
ity measure x 2 . This straightforward fit turns out to 
leave large uncertainties on some nuclear matter proper- 
ties, and with it on nuclear giant resonance excitations. 
Therefore a restricted fit was performed where some bulk 
bulk properties (see last column of table HJ were fixed to 
commonly accepted values [l[ which turn out to yield a 
very good description of giant resonances in 208 Pb. (The 
actual values which were kept fixed can be read off from 
figure [5]) This led to SV-bas where "bas" stands for 
base point because it was used in [ltjj as the origin for a 
systematic variation of bulk properties. 

The new parameterizations take over the twofold strat- 
egy, free fit "min" and bulk restricted fit "bas" . We con- 
sider first the new density dependent term alone leading 
to the parameterizations RD-min and RD-bas. A more 
flexible density dependence is explored by allowing also 
the traditional density dependent term with fixed a = 1 . 
These are the parameterizations RM-min and RM-bas. 
Comparing the RD forces with the RM forces allows to 
learn more about the impact of isovector freedom in the 
density dependence. 

We also compare with results from established SHF 
parameterizations: SkM* as a widely used old standard 
[14| which for the first time managed to deliver accept- 
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able surface energy and fission barriers; SLy6 which had 
been developed with a bias to neutron rich nuclei and 
neutron matter aiming at astrophysical applications [H| ; 
SkI3 adding the freedom of an isovector spin-orbit force 
to obtain an improved description of isotopic shifts of 
r.m.s. radii in neutron rich Pb isotopes (2^1; BSk4 from 
[30| as representative of the series of fits to all available 
nuclei Q having an effective mass of m* /m = 0.92 which 
is best comparable to the new forces here. Each one 
of these parameterizations uses different sets of fit data, 
bias and constraints. Nonetheless, they yield very similar 
results for nuclear bulk properties, but differ in extrap- 
olations and more detailed observables Note that 
these traditional parameterizations are taken as given. 
It is not possible to determine their extrapolations errors 
in the results discussed later on. 

The detailed parameters for the four new parameteri- 
zations are given in tablc|TI]in appendix|B] The extended 
version RM-min and RM-bas have a negligible isoscalar 
parameter b%. A value of zero could be enforced without 
any significant change in quality and other properties. 
However, the isovector parameter b' 3 is large, accompa- 
nied by a similarly large change in the complementing 
63. The extended functional thus exploits its freedom, as 
expected, to tune isovector properties. The parameters 
ai, w and /3j nv are very similar for all four forces. The typ- 
ical value of Qfinv ~ 15 fm 3 means that the functional has 



a simple pole at p = — l/15fm 



-0.067fm" 



This 



means that the functional is not strictly analytical. But 
this pole is safely outside the range of possible single par- 
ticle transition densities (whose absolute value is p{r) but 
which can take any complex phase [13, Ell) an d thus the 
integrations in the projection schemes stay fully in the 
analytical regime. 



III. RESULTS AND DISCUSSION 
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FIG. 1: Average errors on nuclear bulk properties for a va- 
riety of parameterizations as indicated. The errors are ob- 
tained from averaging over the reference set of SV-min [ToT ] . 
Lower panel: errors on binding energy. Upper panel: errors 
on diffraction radius R, surface thickness a, r.m.s. radius r. 



gain is rather small. From the perspective of quality on 
bulk properties the new density dependence is neither 
harmful nor particularly advantageous. It is an equiv- 
alent alternative having the formal advantage to invoke 
only integer powers of density. It remains to been inves- 
tigated how the new parameterizations perform in other 
respects. 



A. Quality with respect to bulk observables 

Figure Q] shows the average errors for the four basic 
nuclear bulk properties. It is to be noted that the older 
parameterizations (SkM*, SLy6, SkI3, and BSk4) were 
tuned with different data and strategies than used here. 
This explains, e.g., the larger energy error for SkI3 which 
was fitted to relative energy deviation and so has more 
bias to light nuclei, or the larger error on the diffraction 
radius R for BSk4 which was not included in its fitting 
set. Nonetheless, comparing with the traditional param- 
eterizations, the progress in accuracy over the years be- 
comes obvious. The parameterizations RD-bas and RD- 
min using the new density dependence provide a descrip- 
tion which is very much comparable to the SV-bas and 
SV-min which were fitted to the same pool, of data but 
using the traditional density dependence of the Skyrme 
functional. The step to RM-min and RM-bas adds more 
flexibility and necessarily should deliver an improved de- 
scription. A reduction of errors is indeed seen, but the 



B. Nuclear matter properties 

Figure [2] shows result for nuclear matter properties. 
The lowest panel shows the incomprcssibility K — 
9 p 2 dp(E/A) where (E/A)(p) is the bulk binding of sym- 
metric matter per nucleon. The predictions are about 
the same for all forces. This isoscalar property is rather 
well determined by the wide span of nuclear sizes in all 
fit data. The extrapolation uncertainty is dramatically 
much larger for the new density dependence (compare SV 
with RD and RM) . This is due to the larger flexibility of 
the new functional which has two parameters instead of 
one in the old density model. The form of density depen- 
dence can have a large influence on the bulk incomprcss- 
ibility. The new functional would allow to adjust much 
different values for K if there is a need to do so. The 
effective mass shows great variations for the traditional 
parameterizations. The larger data sets of the new ad- 
justments (including BSk4) allow to determine its value 
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FIG. 2: Nuclear matter properties (filled boxes) together with 
their extrapolation uncertainty for the selected parameteriza- 
tions. 
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FIG. 3: (Color online) Binding energy per nucleon in ho- 
mogeneous nuclear matter as function of density Lower: for 
symmetric matter. Upper: for neutron matter. Shown are 
results for the "min" parameterizations as indicated. Those 
for the "bas" parameterizations are very similar. The values 
for RD-min are given together with errors bars representing 
the uncertainty in the extrapolation. 



in a rather narrow window close to to* /m = 0.9. The 
"bas" parameterizations do not show error bars because 
the value was fitted. The free fits in the "min" param- 
eterizations confirm the choice with a rather small un- 
certainty. This m* /to = 0.9 is fortunately also the value 
which allows a good description of the giant quadrupolc 
resonance. 

All isovector properties (upper three panels) have large 
uncertainties in the unconstrained fits ("min" series). 
Moreover, their mean value differs visibly from the con- 
strained values which explains the degrading of x 2 when 
switching to the constrained fits (see figure [TJ. It is note- 
worthy that the extended density dependence in RM-min 
leads also to some enhancement of uncertainties in the 
isovector properties. The effect, however, is moderate in 
view of the freedom offered by the 63, b' 3 , terms in the 
density-dependent contribution (|lcjl . 

Figure [3] shows the binding energy curves in homoge- 
neous nuclear matter for SV-min using the traditional 
density dependence and for RD-min as well as RM-min 
with the new functional. All three parameterization pre- 
dict practically the same energies for low and normal nu- 
clear densities, up to about p = 0.28 fm -3 for symmetric 
matter and up to about p = 0.15 fm -3 for neutron mat- 
ter. These are the densities which are explored in finite 
nuclei and it is gratifying to see that the same informa- 
tion entering different models produce the same prcdic- 
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tions in the relevant range of density. However, the ex- 
trapolation to high densities shows significant differences. 
Very large effects are seen for RM-min in neutron mat- 
ter due to the more flexible isovector parameterization. 
The difference in predictive power is nicely signaled by 
the extrapolation uncertainties deduced from the least- 
squares techniques. The error bands remain small in the 
region of low and normal nuclear density and quickly en- 
hance for higher densities. The figure shows only the 
uncertainties for RD-min. The error band has the same 
trend for all parameterizations. It is in magnitude much 
the same for RM-min, but a factor of five smaller for 
SV-min. This difference indicates that the new density 
dependence, having two free parameters instead of one, 
allows more flexibility to adjust the equation of state in a 
broad density range. Unfortunately, the presently avail- 
able empirical data do not allow to exploit this flexibil- 
ity. It may be interesting to invoke results from ab-initio 
calculations to determine the density dependence more 
uniquely. 



Giant resonances 



Figure @] shows the peak frequencies for the three basic 
giant resonances in 208 Pb. The resonance spectra were 
computed with the random-phase approximation done 
self-consistcntly with the same Skyrme interaction as was 
used for the ground state, for technical details see (34H35| . 
The peak positions were deduced from an average over 
the resonance region covering ±2 MeV around the peak. 
The GMR comes out very similar and similarly correct 
for all parameterizations. New are the huge uncertainties 
for all new parameterizations with the rational density 
dependence (RD and RM forces). This results is much 
similar to the results for the incompressibility in figure 
[21 The similarity is not surprising because it is known 
that incompressibility and GMR are closely related to 
each other [3g|. The GDR behaves very similar with old 
and new density-dependence. Unconstrained fits (right- 
most three entries in the middle panel) yield too low peak 
frequencies and large uncertainties while constrained fits 
perform well by construction. There is an intimate rela- 
tionship between GDR and the symmetry energy a sym as 
well as the sum rule enhancement factor k. The larger 
symmetry energies together with lower k in fi.gure[2]corre- 
late with the lower GDR energies in figure [4] The GDR 
in 16 0, however, remains troublesome as it was for all 
previous SHF parameterizations [l(| ED, H3|- The res- 
onance energy stays far below the experimental values. 
Even the extended density dependence using &3,& 3 with 
their larger isovector freedom does not resolve the puz- 
zle. The GQR (upper panel in figure^]) is known to have 
clear relation to the effective mass m*/m (23|. This is 
also found in the present results (when comparing with 
m* jm in figure [5]). The results for the RD and RM forces 
are the same as for the SV forces. This shows that there 
is no influence from density dependence on the GQR. 
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FIG. 4: Average peak energies of giant resonance excita- 
tions in 208 Pb. Lower panel: isoscalar giant monopole reso- 
nance (GMR); middle panel: isovector giant dipole resonance 
(GDR), upper panel: isoscalar giant quadrupole resonance 
(GQR). For the newly fitted forces, we show as error bars 
also the uncertainties in the prediction. The faint horizon- 
tal lines indicate the experimental resonance positions (taken 
from [3IH31). 



0.9 delivers a good reproduction of 



The value m*/m 
the GQR. 



D. Super- heavy elements (SHE) 

Figure[5]shows results for super-heavy elements (SHE). 
The results of the new parameterizations stay close to 
those of SV-bas and SV-min. The new density depen- 
dence makes little difference, mainly a slight growth of 
the uncertainty. Even then the uncertainties remain 
small. As a consequence, the mismatch of binding en- 
ergy in the known SHE 264 Hs persists [H]. We deduce 
from the result that this mismatch is probably not related 
to an insufficient density dependence and thus cannot be 
cured by more flexible density profiles. 

Another important observable in SHE arc fission bar- 
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FIG. 5: Properties of super-heavy elements (boxes) together 
with their extrapolation uncertainty for the selected param- 
eterizations. Lower: binding energy of the nucleus Z = 120, 
N = 180. Middle: Q a value for the nucleus Z = 120, 
TV = 180. Upper: binding energy of 264 Hs (including a cor- 
rection for the rotational zero-point energy [lCj]), the experi- 
mental value is indicated by a dashed horizontal line [381 ]. 



ricrs. Figure [5] shows results for barriers of two groups 
of SHE, at the lower end Rf, Sg, Hs and more heavily 
Z=112,114 all calculated with the new density depen- 
dences, with the traditional one, and compared with ex- 
perimental values [39|,[4(|. The predictions from the three 
forces shown in the figure are very much the same. The 
parameterizations SV-min, RD-min, and RM-min would 
also coincide and are not shown to keep the plot oversee- 
able. The modified density dependence docs not make 
any difference on fission barrriers. The agreement is gen- 
erally satisfying. Note, however, that the experimentally 
observed trend of fission barriers cannot be reproduced 
correctly by any SHF parameterization. The barriers are 
somewhat overestimated at the light side and underesti- 
mated in the heavier group. The same mismatch in the 
trend was found also for all conventional forces like SkI3 
or SLy6 and for fission life-times [12]. The error bars 
indicated in the figure are as large as the deviation, for 
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FIG. 6: (Color online) Fission barriers (filled boxes) of two 
groups of SHE for the selected parameterizations together 
with the experimental values. The extrapolation errors are 
very similar for all three forces. They are indicated for the 
lower group near the left margin and for the upper group near 
the right margin. 



the lighter SHE even larger. This nourishes some hope 
that the trend could be resolved by appropriate retuning 
of the forces. 



E. Energy systematics 

Figure [7] shows the systematics of errors of binding en- 
ergies E^^^ — Eg* p ^ for all known even-even nuclei and 
for the the forces SV-bas, RD-bas, and RM-bas. The re- 
sults for SV-min, RD-min, and RM-min look very similar 
and are thus not shown here. The pattern in figure [7J arc 
practically the same for all three forces: there is perfect 
agreement for the fit nuclei, there is good agreement for 
all nuclei with A < 180 when noting that the soft nuclei 
(the remaining group, neither fit nuclei nor deformed) ac- 
quire additional binding from correlation effects. A sys- 
tematic trend to under-binding of heavy deformed nuclei 
is observed. The mismatch already seen for the case of 
264 Hs in figure [5] was only one spotlight out of this trend. 
Similar as in [lCf we have tried to improve the perfor- 
mance for energies by including deformed nuclei into the 
fit data set. This did improve the situation a bit, how- 
ever, at the price of a much degraded overall quality. We 
find again that density dependence is not the key to the 
solution of this puzzle. 



IV. CONCLUSIONS 

In this paper, we have proposed and investigated an 
alternative form for the density dependent term in the 
Skyrme functional. The traditional form oc p 2+a is re- 
placed by a rational function cx (a- 1ILV p + /3i nv p 2 )/(l + 
a- mv p). The advantage of the new form is that it em- 



8 



RM-bas 



any 

well deformed 
fit nuclei 



RD-bas 




any 

well deformed 
fit nuclei 



SV-bas 




any 

well deformed 
fit nuclei 



50 100 150 200 
mass number A 



250 



ant resonance, super-heavy elements, fission, and energy 
systematics. The predictions in all these observables are 
very similar for both new functionals as compared to the 
recent fit within the standard Skyrme functional. 

Slight differences appear for the uncertainties in the ex- 
trapolations. The new paramcterizations, leaving more 
degrees of freedom, tend to produce larger extrapolation 
errors in some observables. In particular, the errors on 
incomprcssibility and giant monopolc resonance are sig- 
nificantly larger. The extended new parameterization, 
having more isovector flexibility, produces slightly larger 
uncertainties in isovector observables as, e.g., the sym- 
metry energy. 

Large differences shows up for the equation-of-state of 
homogeneous matter at large densities, particularly for 
neutron matter. This is not so surprising because the 
form of density dependence is different. The comple- 
ment of the observation is, in fact, very satisfying. There 
is practically no difference seen in the range of low and 
normal nuclear densities. This means that the adjust- 
ment to finite nuclei has managed to accommodate the 
region of normal densities in spite of the different form 
of the functionals. 

Acknowledgment: This work was supported by BMBF 
under contract no. 06 ER i42D. 



APPENDIX A: DENSITIES AND CURRENTS 



FIG. 7: (Color online) Systematics of binding energies over 
all known even-even nuclei for three parameterizations as in- 
dicated. Fit nuclei are indicated by filled black circles, well 
deformed one by filled red circles, and all remaining one by 
open green circles. 



ploys only integer powers of p and that it is analytical, 
except for one simple pole at p = —ol^ . The new form 
is thus applicable for projection schemes as, e.g., angular 
momentum, center-of-mass projection and particle num- 
ber projection. We have also explored an extended form 
of the functional with one more density dependent term 
allowing for more flexibility in tuning isovector density- 
dependent terms. We have optimized the new function- 
als by adjustment to phenomenological data in standard 
manner and we have explored the properties of the newly 
developed functional in comparison to a couple of tradi- 
tional Skyrme functionals. The straightforward methods 
of least-squares fitting allow also to deduce uncertainties 
on the predictions drawn from the fitted functionals. 

The new forces provide at once an overall quality on 
standard nuclear bulk properties (energy, r.m.s. radius, 
diffraction radius, surface thickness) which is at least as 
good as the one obtained with the standard functional 
using exactly the same fitting data. The extended pa- 
rameterization even provides a slight improvement. We 
have then checked the predictions for a couple of other 
nuclear observables, nuclear matter at equilibrium, gi- 



The densities and currents are defined as 

r q (r) = ]T/ Q ^|V^(r)| 2 

j q (r) = %ml^2f a vl(p+(r)V(p a (r)\ 

Ue? ^ J (AI) 

a£q 
a£q 

£<z(f) = ^2 f a u a v a \ip a (r)\ 2 

a£q 

where q labels the nucleon species with q = p for protons 
and q — n for neutrons. A density without an isospin 
label corresponds to a total density, i.e. p = p p + p n 
and similarly for the other densities. The v a and u a 
are the standard BCS amplitudes, v a for occupation and 
u a the complementing one for non-occupation such that 
v a + U « = 1- The f a is a phase-space weight which serves 
to provide a smooth cutoff of the space of single-particle 
states included in pairing. We use here a soft cutoff 
profile such as, f a = [1 + cxp ((e a - (e F + e cut ))/Ae)]_ 1 
where typically e cut = 5MeV and Ae = e cut /10 (4l1. |42|. 
This works very well for all stable and moderately exotic 
nuclei. For better extrapolation ability away from the 
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RD-bas 


RD-min 


RM-bas 


RM-min 


bo 


-1116.30 


-1116.34 


-1117.47 


-1116.13 


b'o 


-1127.16 


-977.786 


-1911.40 


-1705.29 


h 


15.0000 


14.4966 


15.0000 


10.5256 


b[ 


72.8588 


-18.9894 


71.7330 


16.5372 


62 


110.515 


112.249 


109.873 


109.771 


b' 2 


-284.614 


-317.996 


-81.9695 


-107.028 


h 






-9.64712 


9.95249 








-4997.66 


-7137.87 








1.00000 


1.00000 


h 


4364.93 


4364.74 


4395.28 


4394.86 


b' 3 


3325.54 


4218.63 


15058.7 


16756.7 


Qiinv 


15.1299 


15.1269 


15.2122 


15.3502 


Anv 


16.9485 


17.0715 


16.5759 


17.1909 


64 


63.0960 


55.6340 


64.4335 


61.3840 


64 


36.3802 


52.0100 


33.5588 


38.4206 


Vpair,p 


658.539 


685.826 


708.671 


820.390 


Vpair,n 


601.751 


654.283 


642.556 


740.647 


pO.pair 


0.203601 


0.194093 


0.194622 


0.176469 


2m p 




20.7498207 




h 

2m n 




20.7212601 




e 2 




1.44 




Vtls 




1 






J?Cox 












TABLE II: The parameters of the Skyrme functional (JTc} for 
the four new parameterizations. 

valley of stability, the fixed margin e cut may be modified 
to use a band of fixed particle number cx TV 2 / 3 instead of 



a fixed energy band [43| . 



APPENDIX B: DETAILS OF THE PARAMETERS 

Table |TT] provides the parameters of the four new 
parameterizations discussed in this paper. Dimen- 
sion of length are given in fin and dimensions of 
energy in McV. This means: bo,b' Q <H> [McVfm 3 ], 
bi,b' 1 ,b 2 ,b' 2 o [McVfm 5 ], 63,63,63,63 O [McVfm 6 ], 
64,64 [MeVfm 5 ], a,a inv ,Anv <-> far 3 , e 2 <->• [MeVfm], 
■Sj— <H> [MeVfm 2 ], and rftis is a dimensionless switch fac- 
tor. 

The parameters 6^ and b\ used in the definition (|lc[) 
of the Skyrme functional are chosen to give a most com- 
pact formulation of the functional. The traditional form 
of the Skyrme functional is derived from an effective two- 
body interaction using force parameters U and exchange 
parameters Xi 0, 0, [20, [23| . There is a one-to-one corre- 
spondence of the parameters which reads 



60 


= t (l + \x ) , 




b'o 


= t (§ + Xq) , 




61 


= l[ tl (l + l xl ) 


■f *a(l + \x 2 ) 


b[ 


= i[h$ + xi)- 


t 2 (^ + x 2 )] 


62 


= i [3t x (l + \x x ) 


-t 2 (l+\x 2 


b'2 


= |[3*i(i + ^) 


+■ Hk + x 2 )} 


63 


= 1*3(1 + 1x3) , 




b' 3 


= l M\ + x 3 ) , 
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